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摘 要 考虑 熔 滴 和 小 孔 对 熔 池 的 影响 , 建立 了 基于 FLUENT 软件 的 激光 + 熔化 极 电弧 (GMAW) 复 合 热 源 焊 三 维 瞬 态 熔 
流体 流动 数值 分 析 模 型 . 利用 双 椭 球体 热源 描述 电弧 热 输 入 , 将 激光 热 输入 视 为 热流 峰值 可 调节 的 双 曲 线 旋转 体 热 源 , 其 
热源 分 布 参数 通过 简化 的 小 孔 形 状 尺寸 模型 确定 ; 将 熔 滴 过 渡 过 程 视 为 从 燃 池 上 部 特定 区 域 流 入 熔 池 高 温 液 态 金属 的 过 
程 , 并 通过 建立 液态 金属 流速 对 时 间 的 周期 函数 表征 熔 滴 过 渡 频 率 ; 将 小 孔 视 为 由 激光 致 蒸汽 反作用 力 引 起 的 熔 池 表面 变 
乡 , 以 简化 计算 过 程 , 重点 考虑 小 孔 的 存在 对 熔 池 流体 流 态 的 主要 影响 . 利用 所 建 模型 对 不 同 焊接 条 件 下 的 激光 +GMAW 
复合 热源 焊 小 孔 形 态 、 熔 池 流 体 流 动 和 温度 场 进行 模拟 计算 , 分 析 了 激光 +GMAW 复合 热源 焊 流 场 特征 , 探讨 了 激光 功率 
对 复合 焊 熔 池 动 态 行为 的 影响 规律 . 结果 表明 , 在 1 m/min 焊 速 条 件 下 , GMAW 焊 ( 激 光 功 率 为 0 W) 出 现 驼峰 缺陷 ; 当 激 光 
功率 为 500 W 时 , 驼峰 缺陷 消失 , 但 燃 池 中 无 小 孔 产 生 , 且 流 体 流动 模式 与 GMAW 焊 相 近 ; 而 当 激 光 功 率 增 至 2000 W, 熔 ; 
中 出 现 小 孔 , 使 得 流体 流动 模式 更 为 复杂 . 将 焊 颖 横断 面 形状 尺寸 的 计算 结果 与 实验 结果 进行 比较 ,2 者 吻合 较 好 , 从 而 证 
明了 模型 的 准确 性 和 适用 性 . 
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ABSTRACT Laser+gas metal arc welding(GMAW) hybrid welding fully combines the merits of both laser 
welding and GMAW, which can achieve high quality, high efficiency and comparatively low-cost welding of thin 
and thick plate, thus having great application prospect in manufacturing industry. However, compared with single 
heat source welding, hybrid welding involves more welding parameters and more complicated physical process, 
leading to difficult process optimization. When mismatching the process parameters, welding defects can still ap- 
pear in high-speed welding, which affects the reliability of hybrid welding. Therefore, it is necessary to study the 
physical mechanism in hybrid welding deeply for suppressing welding defects and improving welding stability. In 


hybrid welding, fluid flow in weld pool has a critical influence on the weld formation. So, modeling and simulating 
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激光 + 熔化 极 电弧 (GMAW) 复 合 热 源 焊 综合 了 
激光 焊 和 GMAW 焊 的 双重 优点 , 可 同时 实现 薄板 和 
厚 板 的 高 速 、 优 质 和 相对 低 成 本 的 焊接 , 在 工业 制 
造 领域 具有 巨大 的 应 用 前 景 ". 但 相 较 于 单 热源 焊 。” 析 , 而 汪 任 赁 “同样 采用 二 维 数值 分 析 模 型 对 固定 
接 工艺 , 复合 焊 涉 及 工艺 参数 较 多 , 难于 优化 "1. 若 
工艺 参数 匹配 不 当 , 高 速 焊接 过 程 中 仍 会 出 现 焊接 
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the fluid flow is helpful for understanding the process mechanism completely. To date, however there is only little 
study on Velocity field in hybrid welding due to its complexity. In this work, with considering the effects of droplet 
and keyhole on weld pool, a three dimensional transient model is developed to numerically analyze fluid flow in 
weld pool of lasertGMAW hybrid welding based on FLUENT software. Arc heat input is modeled using an double- 
ellipsoid heat source; laser heat input is regarded as a hyperbolic curve-rotated heat source with changing peak 
power density, its distribution parameters being determined based on the simplified model for keyhole geometry 
and size. Droplet transfer is described as the process of high temperature liquid metal flowing into weld pool from 
the certain domain above the weld pool. Using the built model, the keyhole behavior, fluid flow and temperature 
distribution in lasertGMAW hybrid welding under different welding conditions are calculated. The features of ve- 
locity field in hybrid welding are analyzed and the effect of laser power on the weld pool dynamic behavior is dis- 
cussed. The results show that, in the case of 1 m/min, weld bead hump is generated in single GMAW (laser power 
0 W); when laser power is 500 W, bead hump disappears in welding, but there is no keyhole emerging in hybrid 
weld pool and fluid flow pattern is close to that in GMAW. When increasing laser power to 2000 W, keyhole is 
formed, which makes the fluid flow in weld pool more complicated. The predicted weld geometries and dimen- 
sions for varied laser powers are compared with the measured data, which are in good agreement, thereby indicat- 


ing accuracy and applicability of the established model. 
KEY WORDS hybrid welding, weld pool, fluid flow, keyhole, numerical simulation 
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合 焊 时 熔 滴 进入 小 孔 后 的 扩散 过 程 ; 


= 


缺陷 , 从 而 导致 该 工艺 可 靠 性 相对 较 差 , 阻碍 其 进 
一 步 推 广 应 用 59, 而 焊接 过 程 中 , 熔 池 内 流体 流动 


第 51 卷 


虑 激光 与 电弧 的 相互 作用 , 且 计 算 结 果 缺 乏 实验 验 
正 . Zhou 和 Tsaic 通 过 建立 二 维 数值 模型 , 对 固定 复 
行 了 数值 分 


复合 焊 熔 滴 过 渡 、 电 弧 及 熔 池 进行 了 模拟 计算 .但 2 
者 的 研究 均 没 有 涉及 焊 颖 成 形 , 且 模 型 不 适 于 移动 
热源 焊接 和 三 维 流体 流动 分 析 , 从 而 阻碍 了 模型 的 
进一步 推广 应 用 . 高 志 国 外 对 铝 合 金 复 合 焊 流体 流 


模式 与 驼峰 、 咬 边 及 气孔 等 缺陷 密切 相关 , 对 焊 缝 。、” 动 进行 了 数值 计算 ,但 计算 时 间 较 短 (小 于 0.2 s), 研 
成 形 


热 


i 


源 焊 熔 池 流体 流动 开展 数值 分 析 , 有 助 于 全 面 理 
解 复合 香 焊 接 缺 陷 的 产生 机 理 和 工艺 参数 优化 , 提 


具有 关键 影响 “9. 因此 , 对 激光 +GMAW 复合 


对 焊 缝 成 形 影 响 . 


完结 果 无 法 全 面 地 反映 复合 焊 流 体 流动 特征 及 其 


本 工作 基于 FLUENT 软件 , 建立 了 适用 的 激 


高 焊接 过 程 的 可 靠 性 , 具有 重要 的 学 术 价值 和 工程 ” 光 +GMAW 复 合 热源 焊 三 维 瞬 态 炊 池 流 体 流动 数值 


实用 价值 . 
复合 焊 时 , 激光 致 小 孔 ` 熔 滴 和 烽 池 3 者 之 间 存 


分 析 模 型 , 对 复合 焊 流 体 流 动 模式 进行 了 模拟 计 
算 , 分 析 了 其 流体 流动 特征 , 研究 了 激光 功率 对 复 


在 着 相互 耦合 作用 , 其 熔 池 内 流体 流动 模式 较 单 热 


合 焊 流体 流动 及 焊 颖 成 形 的 影响 规律 , 从 而 为 深入 


源 焊 接 更 为 复杂 中 因此 , 目前 针对 激光 +GMAW 复 
合 热 源 焊 机 理 的 数值 模拟 研究 主要 集中 于 热 场 *" 
及 应 力 场 方面 上 2 和, 对 其 熔 池 流体 流动 研究 非常 缺 
乏 , 且 主 要 集中 于 单一 焊接 条 件 某 一 特定 物理 现象 
的 研究 . Cho 和 Farsong 分 别 采 用 双 李 球体 热源 和 
Gauss 平面 热源 描述 电弧 和 激光 热 输 入 , 基于 
FLOW-3D 软件 对 激光 +GMAW 复合 焊 流 体 流 动 及 
焊 缝 成 形 进 行 了 模拟 计算 , 但 该 研究 仅 适 于 无 小 孔 
产生 条 件 下 的 复合 焊 过 程 . Cho 等 "3 通过 将 Gauss 
平面 热源 与 小 孔 热 源 相 结合 , 建立 了 激光 +GMAW 
复合 热源 焊 三 维 数值 分 析 模 型 , 对 单一 条 件 下 复合 
热源 焊 流 体 流动 模式 进行 了 分 析 , 但 该 模型 没有 考 


技术 支持 和 理论 基础 . 
1 实验 条 件 


理解 激光 +GMAW 复合 焊 内 部 机 理 及 工程 应 用 提供 


实验 母 材 为 304 不 锈 钢板 , 板 厚 为 6 mm; 焊丝 
为 304 不 锈 钢 焊 丝 , 焊丝 直径 为 1.2 mm. 图 1 为 激 


光 +GMAW 复合 热源 焊 示 意图 . 复合 焊 过 程 中 , 激光 


前 置 , 且 与 焊 件 表面 垂直 ; 电弧 在 后 , 与 激光 夹 角 为 
27°, 光 丝 间距 为 3 mm. 焊接 电流 为 170 A, 电弧 电压 


为 20 V, 焊接 速率 分 别 


取 1 和 2 m/min, 激光 功率 分 


别 取 0( 即 为 单 GMAW 焊 ), 500 和 2000 W 3 组 , 离 焦 


量 为 -1 mm, 激光 焦点 直径 为 0.2 mm, 保护 气体 为 


Art+10%CO,( 体 积分 数 ). 
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GMAW wd | si beam 


Laser-wire distance 
1 激光 + 熔化 极 电弧 (GMAW) 复 合 热源 焊 示 意图 
Fig.1 Schematic of laser+gas metal arc welding (GMAW) 


hybrid welding 


2 复合 焊 流 体 流动 数值 分 析 模 型 
2.1 复合 焊 热 源 模型 

2.1.1 电弧 热 输入 - 双 椭 球体 热源 模型 复合 
焊 过 程 中 , 在 电弧 压力 和 激光 致 蒸汽 反作用 力 下 ， 
熔 池 表面 (小 孔 周围 区 域 ) 产 生 一 定 的 变形 , 又 因 复 
合 焊 焊 接 速 度 较 快 , 故 采 用 双 椭 球体 热源 模型 描述 
电弧 热 输入 分 布 模式 . 计算 时 , 热源 中 心 位 于 焊 件 
表面 , 其 热流 密度 分 布 函数 gtkx, 思 习 如 下 所 示 : 


6J3(CAmAD) 
人 0 aibic mT 
w= 37 
串 - ( 2 < 。 -学 (t=0) (1) 
a 6V3(fn0) 
0 abic TNT 


Ch 


| 到 9 3 | (<0) 2 

Q 外 

fitf.=2 (G3) 

式 中 ,x, yz 为 坐标 ; ae as bs, cs 为 热源 分 布 参数 ; 大 和 
大 分 别 为 热源 前 后 2 部 分 热量 的 分 布 系数 ; 为 电弧 
热效率 ; 7 为 焊接 电流 ; U 为 电弧 电压 ; w 为 焊接 速 
率 ; 1 为 时 间 . 

2.1.2 激光 热 输入 -基于 简易 小 孔 模 型 的 锥 体 执 
源 模型 ” 为 了 提高 计算 效率 ,将 研究 重点 放 在 复 
合 焊 流体 流动 数值 分 析 上 , 采用 基于 小 孔 尺 寸 且 热 
流 峰 值 可 调 的 双 曲 线 旋 转 体 热 源 模型 * 描 述 激光 热 
流 分 布 模式 , 其 热流 密度 分 布 函数 gp z) 如 下 所 示 : 

3n.PL 2 
(1 -ee )(4 +B) 


(1 ~X)zz.1 硬 | 37? | (4) 
1 e i e [oa)) 


qi(7,3)= 


Z. 一 各 4 高: 一 就 


ro(z) 二 二 (5) 
Zi—Z 
到 TT 
r.—r 
be ee (6) 
0b 
r.—r 


式 中 , hz 为 坐标 ,r= (v0t) +y ; 1 为 激光 热 效 
率 , 已 为 激光 功率 , ro(z) 为 热源 作用 半径 , x. 和 分别 
为 热源 上 、 下 表面 半径 , z. 和 zz 分 别 为 热源 上 、 下 表面 
z 坐 标 , x 为 热源 上 、 下 表面 热流 峰值 比例 系数 ,a 和 4b 
为 计算 参数 ,4 和 B 为 过 程 计 算 参 数 , 其 具体 计算 公 
式 见 文献 [14]. 
该 模型 综合 考虑 激光 热 输 入 作用 区 域 及 热流 
峰值 沿 焊 件 厚度 方向 的 变化 , 其 作用 半径 沿 焊 件 厚 
度 方向 呈 双 曲线 递减 , 同时 热源 中 心 轴 热流 密度 峰 
值 沿 焊 件 厚度 也 呈 双 曲线 变化 . 当 比 例 系数 x 大 于 1 
时 , 热流 峰值 沿 焊 件 厚 度 方向 按 双 曲 线 模 式 递增 ; 
反之 , 则 递减 . 计算 过 程 中 , 利用 基于 线 热源 的 小 孔 
形状 数值 分 析 模 型 确定 小 孔 尺 寸 , 再 依据 小 孔 尺 寸 
确定 热源 模型 上 、 下 表面 半径 和 高 度 , 其 具体 计算 
步骤 见 文献 [15]; 通过 适当 调节 体积 热源 上 、 下 表面 
热流 峰值 比例 系数 x 间接 考虑 激光 在 小 孔 内 的 多 次 
反射 及 孔 壁 对 激光 能 量 的 Fresnel 吸收 对 激光 能 量 
分 布 的 影响 . 
2.2 熔 滴 过 渡 模 型 

为 了 减少 计算 成 本 , 同时 考虑 熔 滴 对 复合 焊 流 
体 流动 .小 孔 动 态 行为 及 复合 焊 热量 分 布 的 影响 ， 
建立 了 基于 FLUENT 软 件 的 简易 炊 滴 过 渡 模 型 . 计 
算 中 , 将 熔 滴 过 渡 过 程 视 为 液态 金属 以 一 定 速度 由 
熔 池 上 部 焊丝 存在 区 域 流 入 熔 池 的 过 程 , 该 区 域 面 
积 为 焊丝 横 截 面 面 积 , 即 假定 熔 滴 直径 为 焊丝 直径 . 

如 图 2 所 示 , 将 焊 件 求解 区 域 上 部 边界 设 为 液 
态 填 充 金属 的 速度 入 口 (velocity inlet). 在 焊丝 作用 
区 域 (其 面积 即 为 焊丝 横断 面 面 积 ) 设 置 液态 金属 流 
入 的 初始 流速 ( 沿 z 轴 负 方 向 ), 而 其 余 位 置 速度 为 
零 , 即 其 余 区 域 无 液态 金属 流入 , 以 此 表征 焊丝 熔 


Velocity inlet 


Pressure outlet 


Symmetry planes 
X Wall 
图 2 求解 区 域 的 几何 模型 


Fig.2 Geometric model of calculation domain 
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化 后 的 熔 滴 过 渡 . 与 热源 模型 类 似 , 计算 时 焊丝 作 p= en -TY | (10) 
用 区 域 依据 焊接 速度 沿 焊接 方向 移动 . 本 7, 
液态 金属 在 流入 熔 池 过 程 中 仅 受 重力 作用 , 通 大 mL, 1) 
过 适当 提高 液态 金属 初始 的 流入 速度 间接 考虑 电磁 Nh, 


\ 等 离子 流 力 对 炊 滴 过 渡 速 度 的 影响 . 此 外 , 通过 
对 液态 金属 流入 流速 施加 时 间 周 期 函数 , 模拟 燃 滴 
过 渡 频 率 . 假定 在 一 个 熔 滴 过 渡 周 期 内 ,一半 时 间 液 
态 金属 流速 则 为 零 . 计算 时 , 液态 金属 的 初始 流速 设 
定 为 1 m/min, 高 于 文献 [21]; 炊 滴 过 渡 频 率 则 依据 送 
丝 速度 确定 , 而 液态 金属 温度 假定 为 2400 K®. 

2.3 控制 方程 与 边界 条 件 

激光 +GMAW 复合 热源 焊 流 体 流动 非常 复杂 ， 
计算 时 必须 进行 简化 , 故 提出 如 下 假设 : (1) 除 热 导 
率 、 比 热 和 动力 粘度 系数 外 , 材料 其 他 物性 参数 与 
温度 无 关 ; (2) 熔 池 内 流体 为 不 可 压缩 Newton 体 , 流 
动 为 层 流 ; (3) 小 孔 内 压强 为 大 气压 强 . 熔 池 内 流体 
在 电弧 压力 、 表 面 张力 、 电 磁力 、 重 力 、 浮 力 、Maran- 
goni 力 、 激 光 致 蒸汽 反作用 力 等 作用 下 产生 流动 . 描 
述 熔 滴 - 熔 池 -小 孔 传 热 和 流体 流动 的 控制 方程 包括 
能 量 、 动 量 及 质量 守恒 方程 , 其 具体 表达 式 见 文献 
[26], 采用 烩 - 孔 隐 度 法 局 模拟 凝固 燃 化 问题 . 

复合 焊 过 程 中 , 小 孔 对 熔 池 内 流体 具有 强烈 的 
觉 拌 作用 , 数值 模拟 过 程 中 必须 考虑 小 孔 动 态 行为 
对 复合 焊 流 体 流 动 模式 的 影响 . 将 小 孔 视 为 由 激光 
致 蒸汽 反作用 力 、 电 弧 压力 、 表 面 张力 产生 综合 作 
用 下 的 燃 池 表面 变形 , 忽略 小 孔 内 部 等 离子 体 的 热 
场 和 力 场 影响 . 计算 过 程 中 , 熔 池 自由 表面 存在 如 
下 应 力 边界 条 件 : 

Par4p 二 谍 (7) 

式 中 , P 为 作用 于 炊 池 自由 表面 的 压力 , 忆 为 电弧 压 
力 , 闫 蒸汽 反作用 力 , Ps 为 表面 张力 . 
电弧 压力 采用 双 椭 圆 分 布 模式 "1: 


2 3{x— 2 2 
jE ei _ 站 
2T (mi 十 op] Qi j 
(x 三 0) (8) 
2 可 (元 三 2 2 
P(x,y) Sl 2 Le el 一 人 
TT (a ij 十 Op] Qi 
(x<0) (9) 


式 中 ， 上 为 磁 导 率 ; C 为 计算 系数 ; ai ap, Di bi bp 为 
电弧 压力 分 布 参 数 . 

激光 致 燕 汽 反作用 力 是 小 孔 产 生 的 主要 驱动 
力 , 可 表示 为 呈 : 


式 中 , 4。 和 为 计算 系数 , B 为 燕 发 常数 ,及 为 熔 池 
表面 温度 , mm 为 原子 质量 , 广 为 燕 发 潜 热 , N, 为 Avo- 
gadro 常数 , hh 为 Boltzmann 常数 . 

表面 张力 表示 为 : 


Ps=ky (12) 
式 中 , ;为 表面 张力 系数 ,上 为 自由 液 面 曲率 半径 , 其 
计算 公式 见 文献 [26]. 由 于 电弧 压力 、 蒸 汽 反 作用 力 
及 表面 张力 为 表面 力 , 计算 时 需 将 其 转换 为 体积 
力 , 在 熔 池 气 液 界面 处 施加 于 动量 方程 的 源 相 忆 2. 
计算 中 所 需 初 始 条 件 以 及 其 他 能 量 和 动量 边界 条 
件 见 文献 [21]; 而 电磁 力 、 浮 力 等 其 它 作 用 力 按 文献 
[7] 中 的 方式 处 理 . 
2.4 熔 池 自由 表面 

采用 VOF 法 追踪 熔 池 自由 表面 . 流体 体积 函数 
五 的 控制 方程 如 下 式 所 示 : 

ed (13) 

式 中 ,wu, vw 分 别 为 在 x,y,z 方 同 的 流体 速度 . 
2.5 求解 方法 

采用 流体 动力 学 软件 FLUENT 对 复合 焊 流体 
流动 进行 模拟 计算 . 由 于 对 称 , 取 一 半 工 件 作 为 计 
算 区 域 , 尺寸 为 50 mmx20 mmx9 mm, 上 方 3 mm 为 
空气 (图 2). 利用 GAMBIT 软件 建立 计算 区 域 的 几 
何 模型 , 并 对 其 划分 网 格 . 对 于 空气 分 布 区 域 , 除 对 
称 面 外 , 上 部 边界 设 为 速度 入 口 , 其 余 边界 设 为 压 
力 出 口 ; 而 对 于 金属 分 布 区 域 , 除 对 称 面 外 , 其 余 3 
个 边界 设 为 墙壁 . 利用 用 户 自 定义 函数 定义 边界 条 
件 (包括 高 温 液 态 金属 流入 熔 池 初始 速度 ) 材料 属 
性 以 及 能 量 和 动量 的 源 项 . 计算 时 , 综合 考虑 计算 
精度 和 效率 两 方面 因素 , 本 文采 用 变化 的 时 间 步 
长 , 其 取 值 范围 为 104 人 10*s; 并 通过 PISO 算法 ( 求 
解压 力 的 隐 式 分 离 算法 ) 计 算 流 场 . 
3 复合 焊 流体 流动 的 数值 计算 

利用 上 述 模 型 , 通过 FLUENT 软件 对 不 同 焊接 
条 件 下 激光 +GMAW 复合 热源 焊 流 体 流 动 进行 模拟 
计算 . 计算 中 所 采用 的 材料 物性 参数 及 计算 参数 见 
文献 [23,26]. 依据 上 述 方法 确定 模型 分 布 参数 , 对 于 
单 GMAW 焊 ( 激 光 功 率 为 0 W, 焊 枪 仍 与 焊 件 表面 呈 
一 定 角度 ), 模型 分 布 参数 为 : a=4 mm, a=6 mm, b= 
4 mm, c=2 mm, an=3 mm, ap=4 mm, b=3 mm. 其它 


工艺 条 件 下 模型 分 布 参 数 同样 依据 上 述 方法 确定 . 
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图 3 为 单 GMAW 焊 众 池 纵 截面 流体 流动 的 计 
算 结果 . 可 以 看 出 , 由 于 焊接 速度 较 快 , 在 电弧 压力 
作用 下 , 炊 池 前 部 流体 迅速 流向 熔 池 后 部 , 最 大 流 
速达 2.5 m/s, 而 在 熔 池 前 部 电弧 作用 区 域 形 成 较 大 
凹陷 ; 凹陷 前 壁 处 液态 金属 层 非 常 薄 , 该 位 置 液态 
金属 以 较 大 速度 流向 凹陷 底部 . 当时 间 拓 0.3 s 时 ， 
由 于 液态 熔 池 体积 较 小 , 长 度 较 短 , 熔 池 四 陷 下 部 
的 液态 金属 在 电弧 压力 排挤 下 , 由 熔 池 底部 向 后 流 
动 , 在 熔 池 中 部 紧 靠 热源 处 产生 顺 时 针 涡 流 ; 而 在 
炊 池 后 部 , 受 表面 张力 和 流体 静 压 力 影响 , 产生 逆 
时 针 涡流 , 如 图 3a 所 示 . 随 着 电弧 热源 的 移动 , 液态 
金属 逐渐 在 熔 池 后 部 堆积 , 并 在 金属 堆积 区 域 与 炊 


前 部 液态 金属 向 后 部 的 流动 以 及 堆积 区 液态 金属 向 
熔 池 前 部 的 回流 , 使 得 熔 池 前 部 出 现 新 的 金属 堆积 ， 
如 图 3c 所 示 . 而 由 图 3d 可 知 , 当 太 2 s 时 , 由 于 过 渡 区 
域 温度 较 低 , 由 熔 池 前 部 流入 熔 池 后 部 液态 金属 较 
少 ; 在 流体 静 压 力作 用 下 , 熔 池 后 部 液态 金属 分 别 
向 燃 池 后 部 和 前 部 流动 , 熔 池 前 部 则 重新 出 现 逆 时 
针 涡 流 . 当 过 渡 区 域 液态 金属 凝固 后 , 焊 颖 中 便 形 
成 驼峰 缺陷 . 由 于 焊接 时 , 熔 池 内 流体 流动 受 多 种 
因素 影响 , 其 熔 池 流 态 易 产生 波动 , 文中 给 出 的 仅 
为 某 一 特定 时 刻 的 流体 流动 模式 . 图 4 为 该 工艺 条 
件 下 GMAW 焊 焊 缝 上 表面 计算 结果 与 实验 结果 . 可 
见 , 计算 结果 能 够 较 好 地 反映 焊 颖 成形 , 从 而 证 明 


池 前 部 四 陷 之 间 形 成 一 个 注 层 金属 过 渡 区 域 ; 该 区 
域 液 态 金属 高 速 流向 熔 池 后 部 , 并 克服 流体 静 压 力 ， 
在 熔 池 后 部 形成 逆 时 针 涡 流 , 如 图 3b 所 示 . 由 于 过 渡 
区 域 液态 金属 热量 较 少 , 凝固 速度 较 快 ; 而 且 该 区 域 
热 毛 细 稳 定 较 差 , 易 产 生 横 向 收缩 ,从 而 阻碍 炊 池 


了 模型 的 准确 性 . 需要 指出 , 图 3 所 示 流 体 流动 模式 
是 本 研究 条 件 下 产生 驼峰 缺陷 的 关键 因素 , 而 流体 
流速 由 于 受 焊 速 、 焊 接 电流 等 多 种 因素 影响 , 故 其 
峰值 与 驼峰 缺陷 形成 的 关系 不 明显 . 此 外 , 驼峰 缺 
陷 的 形成 还 受 沿 焊 缝 方向 液态 金属 表面 张力 的 稳 
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3 不 同时 刻 GMAW 焊 熔 池 纵 截面 的 温度 场 和 流 场 分 布 
Fig.3 Temperature and flow field distributions at longitudinal section of weld pool in GMAW for 0.3 s (a), 0.6 s (b), 1.1 s 


(c) and 2 s (d) with laser power P.=0 W (The white lines with arrows are used to illustrate the calculated fluid flow 


patterns more clearly) 
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4 GMAW 焊 颖 上 表 


| 


的 实验 和 计算 结果 


Fig.4 Experimental (a) and calculated (b) results of weld bead top surface in GMAW 


定性 、 材 料 性 能 等 因素 影响 中. 

5 为 x=21 mm 处 GMAW 焊 熔 池 横 截面 温度 
场 及 流 场 的 计算 结果 . 可 见 , 当 太 0.5 s, 焊 件 被 电弧 
热量 熔化 , 并 在 电弧 压力 作用 下 产生 熔 池 表面 变 
形 , 液态 金属 流向 熔 池 底部 ; 当 太 0.6s 时 , 熔 池 表面 
凹陷 被 液态 金属 填充 , 对 称 轴 O=0) 处 液态 金属 由 熔 
池 表 面向 下 流动 , 这 与 电磁 力 影响 有 关 ; 而 熔 池 表面 
离 对 称 轴 稍 远 处 流体 在 表面 张力 作用 下 向 熔 池 边界 
流动 . 随 着 热源 的 持续 移动 , 电弧 压力 和 电磁 力 对 该 
位 置 的 影响 逐渐 减 小 ; 当 太 1.2 s, 金属 隆起 形成 , 并 
在 表面 张力 及 浮力 作用 下 , 产生 顺 时 针 涡流 , 如 图 5c 
所 示 . 当时 间 进 一 步 推 移 , 该 位 置 流 体 流速 减 小 , 熔 
池 开 始 凝固 , 形成 焊 缝 余 高 (图 5d). 

6 和 7 分 别 给 出 了 激光 功率 为 500 W 时 复合 
焊 熔 池 温 度 场 和 流 场 的 计算 结果 . 由 于 激光 能 量 密 
度 较 小 , 燕 气 反作用 力 影 响 不 明显 , 复合 焊 熔 池 中 
无 小 孔 形 成 ; 与 单 GMAW 焊 相 似 , 焊接 时 仅 在 熔 池 
前 部 电弧 作用 区 域 产 生 较 大 熔 池 变形 , 而 且 熔 池 纵 
截面 和 横 截 面 流体 的 基本 流动 模式 也 与 GMAW 焊 
相近 . 但 与 单 GMAW 焊 相 比 , 当 激 光 功 率 为 500 W 
时 , 复合 焊 熔 池 中 无 明显 金属 堆积 现象 , 驼峰 缺陷 
消失 , 这 可 能 是 由 于 激光 能 量 的 加 入 使 焊接 热 输 入 
增加 , 液态 熔 池 体积 增 大 , 冷却 时 间 有 所 增长 , 从 而 
为 熔 池 后 部 液态 金属 的 回流 提供 了 较为 充足 的 时 
间 ; 同时 , 焊接 热流 密度 的 增加 使 得 熔 池 内 液态 金 
属 的 温度 升 高 , 动力 粘度 减 小 , 流动 性 增强 , 利于 改 
善 焊 缝 成形 . 此 外 , 前 置 激 光 对 焊 件 的 预 热 , 改善 了 
液态 金属 在 焊 件 表面 润 湿性 , 利于 其 横向 铺展 , 也 
对 驼峰 缺陷 的 产生 具有 一 定 抑制 作用 . 可 以 看 出 ， 
即使 功率 值 较 小 时 , 激光 的 存在 仍 能 较 好 得 抑制 驼 
峰 缺 陷 , 但 对 熔 池 内 流体 流动 的 基本 模式 影响 较 
小 , 仅 改 变 了 热流 分 布 特征 . 

8 给 出 了 激光 功率 为 2000 W 时 复合 焊 燃 池 纵 
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图 5 不 同时 刻 GMAW 焊 熔 池 横 截面 的 温度 场 及 流 场 


分 布 
Fig.S Temperature and flow field distributions of weld 
pool at cross section in GMAW for 0.5 s (a), 0.6 S 
(b), 1.2 s (c) and 1.5 s (d) with P.=0 W and x=21 mm 


此 种 工艺 条 件 , 在 六 0.3s 和 大 1.0s 时 的 小 孔 深 
度 较 万 0.6s 时 大 , 即 复合 焊 过 程 中 小 孔 出 现 振荡 ， 
并 非 稳 定 存 在 , 这 是 因为 复合 焊 液态 熔 池 体积 较 
大 , 且 存 在 熔 滴 过 渡 , 同时 电弧 压力 、 熔 滴 冲 击 力 和 
电磁 力也 对 熔 池 内 流体 流动 产生 较 大 影响 , 因此 ， 
小 孔 壁 尤其 是 小 孔 后 壁 受 力 难 于 平衡 , 从 而 导致 小 
孔 无 法 稳定 存在 , 小 孔 深 度 出 现 振荡 . 但 由 图 8 可 
知 , 其 深度 的 振荡 范围 基本 稳定 . 此 外 , 由 于 复合 焊 
焊 速 较 快 , 小 孔 并 非 关 于 激光 中 心 轴线 对 称 , 而 是 
明显 向 后 方 倾斜 (图 ga 和 bj, 体现 出 了 高 速 焊 过 程 
中 弯曲 小 孔 形 态 的 主要 特征 外. 而 对 于 弯曲 小 孔 , 采 


截面 的 温度 场 和 流 场 . 可 以 看 出 , 由 于 激光 功率 较 
大 , 激光 致 蒸气 反作用 力 影响 增加 , 焊接 时 熔 池 内 出 
现 小 孔 ; 激光 能 量 的 加 入 不 仅 消除 了 焊 颖 的 驼峰 缺 
陷 , 而 且 使 复合 焊 炊 深 明显 增加 , 表明 在 激光 功率 较 
高 时 , 激光 热 输 入 对 复合 焊 熔 深 具 有 关键 影响 . 


上 轴 对 称 的 激光 体积 热源 , 有 利于 将 大 部 分 能 量 稳 
定 的 作用 于 小 孔 前 壁 , 使 其 受 力 稳定 性 相对 较 大 ， 
这 与 高 速 激光 深 焙 焊 实 际 烛 接 物理 过 程 相 符 ”” 
可 见 , 模型 能 够 反映 出 小 孔 动态 形态 的 主要 特征 
体现 其 对 复合 焊 流体 流动 的 主要 影响 
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图 6 激光 功率 为 500 W 时 不 同时 刻 复合 焊 熔 池 纵 截面 的 温度 场 和 流 场 分 布 
Fig.6 Temperature and flow field distributions at longitudinal section of hybrid weld pool for 0.3 s (a), 0.6 s (b), 1.0 s (c) 


Color online 
图 7 激光 功率 为 500 W 时 不 同时 刻 复 合 焊 熔 池 横 截面 温度 场 和 流 场 分 布 
Fig.7 Temperature and flow field distributions at cross section of hybrid weld pool for 0.5 s (a), 1.0 s (b), 1.1 s (c) and 1.3 s 
(d) with P.=500 W and x=21 mm 


进一步 分 析 可 知 , 复合 焊 时 , 在 电弧 压力 作用 
下 , 小 孔 口 附近 熔 池 表面 产生 一 定 的 下 凹 变形 ; 而 
受 蒸 气 反 作用 力 影 响 , 小 孔 壁 附近 流体 流速 较 大 . 
对 于 小 孔 前 壁 , 由 于 其 受到 相对 稳定 的 蒸汽 反作用 
力作 用 , 其 附近 液态 金属 流速 明显 大 于 其 他 位 置 ， 
最 高 流速 可 达 5 m/s 以 上 , 而 且 该 位 置 液 态 金 属 始终 


由 小 孔 口 向 下 流动 , 再 经 小 孔 底部 流向 熔 池 后 方 ， 
然后 在 浮力 作用 下 向 上 流动 ;上 浮 金 属 一 部 分 向 炊 
池 尾 部 流动 , 一 部 分 向 炊 池 前 部 回流 , 并 在 熔 池 中 
部 形成 顺 时 针 涡 流 . 对 于 小 孔 后 壁 , 除 受 蒸汽 反 作 
用 力 外 , 熔 滴 冲击 力 \ 电 弧 压 力 以 及 流体 动 压力 和 
静 压 力 对 其 影响 相对 较 大 , 故 其 形态 稳定 性 较 差 ， 
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8 激光 功率 为 2000 W 时 不 同时 刻 复合 焊 熔 池 纵 截面 的 温度 场 和 流 场 分 布 
Fig.8 Temperature and flow field distributions at longitudinal section of hybrid weld pool for 0.3 s (a), 0.6 s (c), 1.0 s (b) 


and 1.4 s (d) with P.=2000 W 


其 附近 液态 金属 流 态 也 较为 复杂 , 与 小 孔 形态 有 
关 . 当 小 孔 较 深 时 , 小 孔 后 壁 下 部 液态 金属 所 受 蒸 
汽 反 作用 力克 服 其 表面 张力 , 向 下 流动 , 并 汇 入 将 
池 底 部 向 后 流动 的 金属 ; 而 上 部 金属 则 受 电弧 压力 
影响 较 大 , 产生 向 上 的 流动 , 再 经 炊 池 表面 流向 炊 
池 后 部 ; 但 该 部 分 金属 在 离 小 孔 稍 远 处 , 受到 熔 池 
冲击 力 、 电 磁力 以 及 炊 池 后 部 回流 金属 影响 , 又 向 
熔 池 底部 流动 , 然后 再 流向 炊 池 后 部 , 如 图 8a 和 c 所 
示 ; 而 当 小 孔 较 浅 时 , 其 附近 液态 金属 仅 存在 沿 小 
孔 后 壁 向 上 的 流动 , 如 图 8b 所 示 , 这 与 受到 向 下 流 
动 的 小 孔 前 壁 金属 的 排挤 及 电弧 压力 作用 有 关 . 
可 见 , 当 激 光 功 率 为 2000 W 时 , 由 于 小 孔 的 存 
在 , 使 得 复合 焊 纵 截面 熔 池 内 流体 流动 模式 与 单 
GMAW 焊 和 小 功率 激光 复合 焊 有 很 大 不 同 ; 在 此 种 
工艺 条 件 下 , 炊 池 中 存在 相对 稳定 的 顺 时 针 流 动 ， 
这 将 有 利于 把 更 多 电弧 热量 带 入 次 池 底部 , 增 大 次 
深 , 同时 还 有 助 于 填充 金属 中 合金 元 素 在 炊 池 内 的 


均匀 分 布 . 

同样 , 小 孔 的 搅动 作用 使 得 复合 焊 熔 池 横 截面 
形 貌 及 流体 流 态 发 生 了 变化 . 图 9 为 激光 功率 为 
2000 W 时 复合 焊 熔 池 横 截面 流 场 的 计算 结果 . 可 以 
看 出 , 当 太 0.44 s, 受 电弧 电流 影响 , 工件 金属 开始 熔 
化 , 并 在 电弧 压力 作用 下 , 产生 一 定 的 表面 变形 ; 而 
当 友 0.5 s 时 , 在 蒸汽 反作用 力作 用 下 形成 较 深 的 小 
孔 , 熔 池 中 流体 自 上 而 下 流动 , 熔 池 横断 面 形状 具 
有 激光 深 熔 焊 特 点 . 随 着 焊接 热源 的 移动 , 该 位 置 
小 孔 被 液态 金属 填充 , 小 孔 消 失 , 且 燃 池上 部 宽度 
受 电流 热流 影响 而 增加 . 当 大 0.8 s, 余 高 形成 , 熔 池 
中 流体 流速 明显 降低 . 熔 池 上 部 宽度 较 大 , 而 下 部 
较 窗 , 即 熔 池 横断 面 同 时 具有 GMAW 焊 和 激光 深 熔 
焊 特 点 . 此 时 , 受 电 磁力 和 熔 池 冲击 力 影 响 , 对 称 轴 
中 心 处 流体 自 上 而 下 流动 . 而 当 太 1.2 s, 该 位 置 受 电 
弧 压 力 及 电磁 力 影响 较 小 , 在 浮力 作用 下 , 对 称 轴 
处 流体 自 下 而 上 流动 , 并 受 表面 张力 影响 , 在 熔 池 
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成 顺 时 针 涡 流 . 见 , 由 于 焊 速 大 幅 提 高 , 复合 焊 熔 深 明 显 减 小 ; 同 
要 指出 , 激光 对 驼峰 缺陷 的 抑制 作用 受 焊 接 时, 激光 对 熔 池 的 热 影响 相对 减弱 , 焊 缝 中 再 次 出 


速度 影响 . 图 10 给 出 了 激光 功率 为 500 W, 焊接 速 现 一 定 程 度 的 驼峰 缺陷 . 


率 为 2 m/min 时 复合 焊 熔 池 纵 截面 的 计算 结果 . 可 11 和 12 分 别 给 出 了 激光 功率 为 500 和 2000 W 
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Color online 
9 激光 功率 为 2000 W 时 不 同时 刻 复合 焊 炊 池 横 截面 的 温度 场 和 流 场 分 布 
Fig.9 Temperature and flow fields at cross section of weld pool in hybrid welding for 0.44 s (a), 0.5 s (b), 0.8 s (c) and 1.2 s 
(d) with P.=2000 W and x=20 mm 
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10 激光 功率 为 500 W, 焊接 速率 为 2 m/min 时 复合 焊 熔 池 纵 截面 的 温度 场 和 流 场 分 布 
Fig.10 Temperature and flow field distributions at longitudinal section of hybrid weld pool with P.=500 W and welding 


rate 2 m/min 
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11 复合 焊 焊 颖 上 表面 实验 与 计算 结果 的 比较 
Fig.11 Comparison of experimental (a, ¢, e) and calculated (b, d, f) weld bead top surfaces in hybrid welding with P.=500 W 
at welding rates 1 m/min (a, b) and 2 m/min (c, d) and P.=2000 W at 1 m/min (e, 全 
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12 复合 焊 焊 缝 横 截面 实验 和 计算 结果 的 比较 
Fig.12 Comparison of experimental (a, c) and calculated (b, d) weld cross sections in hybrid welding with P.=500 W (a, b) 
and P.=2000 W (e, d) 


时 复合 焊 烛 颖 上 表面 和 横 截面 实验 与 计算 结果 的 。” 陷 消 失 ; 而 当 焊 速 增 大 到 2 m/min 时 , 驼峰 重新 出 现 . 
比较 . 其 中 , 焊 颖 余 高 依据 填充 金属 的 隆起 获得 , 而 (5) 当 激 光 功 率 为 2000 W, 焊 速 为 1 m/min 时， 
焊 颖 横断 面 余 高 及 爆 颖 形状 尺寸 的 具体 计算 方法 ”复合 焊 熔 池 中 出 现 小 孔 , 但 其 深度 不 稳定 , 始终 处 
见 文 献 [16]. 可 以 看 出 , 2 者 吻合 较 好 , 仅 在 熔 合 线 走 。” 于 震荡 过 程 中 ; 小 孔 前 壁 附近 金属 流速 较 大 , 且 始 


向 方面 存在 一 定 误差 , 从 而 证 明了 所 建 模型 的 准确 终 由 小 孔 口 流向 小 孔 底部 ; 小 孔 后 壁 金 属 的 流 态 与 
性 和 适用 性 . 小 孔 动态 行为 有 关 . 小 孔 的 存在 对 熔 池 内 流体 具有 
结论 较 大 搅动 作用 , 熔 池 中 部 始终 存在 顺 时 针 涡流 , 利 
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